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Abstract: In this paper we introduce Redberry — an open source computer algebra system 
designed to manipulate with symbolic tensorial expressions. It implements basic computer alge- 
bra system routines as well as complex tools for real computations in physics. Redberry core 
provides common for majority of computer algebra systems tools for expressions manipulation, 
generalized on tensorial objects, as well as tensor-specific features: indices symmetries, WTf^- 
style input, natural dummy indices handling, multiple index types etc. The high energy physics 
package includes tools for Feynman diagrams calculation: Dirac and SU(N) traces, Levi-Civita 
simplifications and tools for one-loop calculations in general field theory. In the present paper we 
give detailed description of Redberry functionality: from basic manipulations with tensors to real 
Feynman diagrams calculation, accompanied by many examples. We also introduce graph repre- 
sentation of a tensor — the basic underlying idea of the Redberry architecture, which clarifies 
a deep connection between symbolic tensor algebra and computational graph theory. Redberry 
is written in Java 7 and provides convenient Groovy-based user interface and extensive API for 
developers. The source code of Redberry is covered by more then 800 tests. Redberry is available 
from redberry.cc. 
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1 Introduction 



1.1 Background 

General-purpose computer algebra systems (CASs) have become an essential part of many 
scientific calculations. At the moment there is no CAS that provides generic tools for calcula- 
tions both with ordinary symbolic expressions and tensors, and at the same time have general 
programming capabilities. Among the existing standalone systems we would like to emphasize 
Cadabra developed by Kasper Peeters [1, 2]. In our opinion this CAS provides the most com- 
prehensive set of tensor-oriented tools: tensor symmetries, multiple index types, dummy indices 
handling, substitutions, comparison of tensorial expressions and tensor-specific transformations 
like contractions with metrics, symmetrization, simplification of tensorial expressions etc. On 
the other hand, Cadabra does not provide general programming environment, i.e. such basic con- 
structions as loops, branching and low-level expression manipulation. This makes it impossible 
to develop user-defined algorithms and modules in Cadabra^ . Also, there is a number of software 
packages which deal with tensors, built on top of general-purpose CASs and aimed at solving par- 
ticular problems in physics and mathematics. Among these packages we would like to highlight 
FeynCalc [3], which provides facilities for calculations in high energy physics (in particular for 
Feynman graphs). The main drawback of such software packages is that they consider tensors 
as additional but not a fundamental entities, which restricts the scope of their application by 
the particular set of problems for which they are intended. Authors have no purpose to make a 
review of all tensor-oriented software here but only to emphasize the most relevant works. 

The main challenge of Redberry is to develop the computer algebra system with native support 
of tensorial objects, i.e. to generalize basic routines of computer algebra systems onto tensors. 
Such an approach leads to a fundamentally different point of view on the architecture of the 
CAS, because tensorial expressions have principally different structure and manipulation rules 
in contrast to indexless expressions. Redberry is developed in such a way that its functionality 
is available through the modern high-level general-purpose programming language with excellent 
IDE support. Thereby, Redberry gives a common language for both particular calculations and 
development of new modules and algorithms. As an example of the successful use of such a 
model, one can mention SymPy [4] computer algebra system, which uses Python as a main user 
interface. 

The key features of Redberry include: support of permutational symmetries of tensor indices, 
extensive tools for comparison of tensorial expressions, native dummy indices handling (including 
automatic clash resolution), multiple index types, I^TgX-style input/output and a comprehensive 
set of tensor-specific transformations. On top of core functionality, Redberry implements a 
number of tools for calculations in high-energy physics including calculation of Feynman diagrams 
and one-loop counterterms in general field theory. 

Redberry core is written in Java, while the user interface is written in Groovy and is intended 
to be used within the Groovy environment. Redberry does not provide GUI by itself, but modern 
open source IDEs^ give a very convenient way of Redberry usage, including syntax highlighting 
and code completion. Installation instructions and additional documentation can be found on 
the Redberry website: http://redberry.cc. 

1.2 The introductory example 

Let's start from the minimal but still complex real example to quickly dive into Redberry (the 
next section covers every aspect in details). The example demonstrates calculation of differential 
cross section of the Compton scattering in scalar QED: 



^Though it is available on the C++ level of Cadabra. 

^we recommend JetBrains IntelliJIDEA because of its brilliant support of Groovy 
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Listing 1 : Compton scattering in scalar QED 



1 //photon-scalar-scalar vertex 

2 def VI = 'V_i[p_a, q_b] = -Be* (p_i+q_i) ' . t 

3 //photon-photon-scalar-scalar vertex 

4 def V2 = 'V_{ij} = 2*I*e**2*g_{ i j } ' . t 

5 //scalar propagator 

6 def P = 'D[k_a] = -1/ (k''»k_a-m**2) ' .t 

7 //matrix element 

8 def M = {'M^ij =' 

9 + 'Vi[pl_a,pl_a+kl_a>D[pl_a+kl_a]*V''j[-p2_a,-pl_a-kl_a] ' 

10 + ' +Vj [pl_a, pl_a-k2_a}tD [pl_a-k2_a]*V'^i [-pl_a+k2_a, -p2_a] ' 

11 + ' H-V^i j') .t 

12 //substituting vertices and propagator in matrix element 

13 M = (VI & V2 & P) » M 

14 //squared matrix element 

15 //here minus is due to complex conjugation 

16 def M2 = M » 'M2 = -M_ij*M^ij' .t 

17 //expand squared matrix element and 

18 //eliminate metrics and Kronecker deltas 

19 M2 = (ExpandAll & EllminateMetrlc^ >> M2 

20 M2 = 'd_i^i = 4'.t » M2 

21 //defining mass shell and Mandelstam variables 

22 def mandelstamSubstitutions = setMandelstam{ 

23 ['kl_i': '0', 'pl_i': 'm', 'k2^i': '0', 'p2_i': 'm']) 

24 //substituting mass shell and Mandelstam definitions 

25 M2 = mandelstamSubstitutions >> M2 

26 M2 = 'u = 2*m**2-s-t' .t >> M2 

27 //factor terms 

28 M2 = Factor » M2 

29 println M2 



This script will print the well-known expression for the squared matrix element of the Compton 
scattering in scalar QED: 



which becomes equal to the differential cross section da/dCl after multiplying by l/(647r^s). 

The code above speaks for itself, so leaving aside the physical aspects of the problem, let's 
focus only on programming aspects. t construct converts a string representation into 

a computer object. The input notation for tensors is the same as used in IM^^with minor 
syntax modifications^ (see Sec. 2.1). Redberry uses the predefined notation for some built-in 
tensors, such as I for image one, or g_ij and d^i_j for the metric tensor and Kronecker 
delta respectively. Thus, for example, when we perform contractions with metrics and deltas in 
the line 19, Redberry automatically takes into account that Qai = . However, since the 
dimension of the space affects Kronecker trace value and it is not specified anywhere explicitly, 
we substitute the trace manually (on line 20). 

As could be seen from the example, state of index (upper or lower, i.e. contravariant or 
covariant) is important — indices are considered to be contracted if they have different states. 

All transformations in Redberry (e.g. substitutions, Expand, EliminateMetrics etc.) are 
first-class objects and can be assigned to variables (like in line 22). They can be applied to 

^One can omit curly braces where it does not cause ambiguity. 



\M\' 



(s — m?)"^ {—s — t + m?)"^ 



+ 2m^st - Arr?sh + 2s^t + mH'^ + sh^) 
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mathematical expression using » operator. Detailed discussion on Redberry transformations 
usage and list of basic built-in transformations can be found in Sec. 3. 

The last point that should be discussed here is setMandelstam ( . . .) function, which is 
used in the lines 22-23. It takes Lorentz four-vectors with the corresponding particle masses and 
simply returns a list of substitutions following from provided values. This function is a part of 
redberry-physics (see Appx. B). 

1.3 Paper structure 

The whole paper is accompanied by many Groovy examples. Each example can be executed 
by simple wrapping the code in the following way: 



@Grab (group = ' cc . redberry ', module = 'groovy', version = '1.1') 

import cc . redberry . groovy . Redberry 

import static cc . redberry . core . tensor . Tensors . * 
import static cc . redberry . groovy . RedberryPhysics . * 
import static cc . redberry . groovy . RedberryStatic . * 

use (Redberry) { 

//example code 

} 



This will automatically download Redberry and all required dependencies. Some examples also 
may require additional import statements, in this case it will be noted in a footnotes. 

In Sec. 2 we give description of Redberry syntax and basic functionality, needed for usage and 
understanding of main Redberry principals. In Sec. 3 one can find a list of selected transforma- 
tions and general aspects concerning their application. Sec. 4 illustrates Redberry application in 
high energy physics: complete calculation of Compton scattering in QED and complete calcula- 
tion of one-loop counterterms of non-minimal vector field and minimal fourth order operator. In 
Sec. 5 we consider selected architectural solutions used in Redberry including graph representa- 
tion of tensors. 

2 Basics 

2.1 Basic input/output 

Let's start detailed description of Redberry from discussion on syntax used to define tenso- 
rial expressions. As mentioned above, there is a special syntax construction to convert string 
representation of tensor into computer object: ' . . . ' .t ^. Here is a usage example: 

def tensor = 'F'^{A}_{B \\mu \\nu} * a'.t 

The general conventions on the expression input are the same as in many other CASs. The main 
difference arises from the tensorial nature of Redberry. The following code gives an idea of the 
valid Redberry syntax for tensors: 



1 def t 

2 //same, but without braces 



In the Java API this conversion is performed using static method Tensors .parse (String) 
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3 t = 'F*A_B\\mu\\nu * a'.t 

4 //here braces are necessary 

5 t = 'F'^{A}_{B \\mu \\nu} * a'.t 

6 //here braces are necessary 

7 t = 'F^A_{B_{21}C\\mu\\nu} '.t 

8 //tensor name in the LaTeX notation 

9 t = ' \\Gainma^A_{\\alpha\\beta} '.t 

10 //complete LaTeX style 

11 t = ' \\Gainma{ }^A{ }_{\\alpha\\beta} ;t 



As one can see, in many cases curly braces can be omitted when inputting the indices of tensor. 
Braces are necessary when one needs to separate indices by spaces, or input a subscripted index. 
Arbitrary types of indices are not ahowed in the current version of Redberry: there are 8 pre- 
defined types including Latin and Greek indices (both lower and upper cases) and indices with 
strokes. For further information about indices see Sec. 2.5. In contrast, there is no restrictions 
on tensor names, so we can input a tensor with Greek name written in the IM^i^X notation, like 
in line 9 or 11 (both lines will produce the same tensor). 

The notation for functions (both user defined and built-in scalar functions: sin, cos etc.) is 
the same as in Wolfram Mathematical 



//user defined tensor field 

t = 'F^a.bcdCC^a.bc, p'^a]'.t 

//some Redberry built-in functions 

t = 'Sin[m**2 - p_m*p^m] - Log[x/2]'.t 

t = 'Power[a, b] '.t //same as 'a**b'.t 



It is important to note, that according to the meaning of scalar functions their arguments must 
always be scalar, so the following code will exit with error: 

def t = 'Sin[A_m] '.t 

\> IllegalArgumentException 

All Redberry objects have a default string representations. So, in order to output some 
expression to the stdout one can simply use print In keyword in Groovy scripts. Additionally 
to the default Redberry output format, there are several others. The following example illustrates 
their usage^: 



1 deft = 'F_mn'^{\\alpha\\beta}/ (a+b) '.t 

2 //default output format 

3 println t 

[> (b+a) ** (-1) *F_{mn} { \alpha\beta} 

4 //LaTeX format 

5 println t . toString (LaTeX) 

[> \frac{ 1 } { (a+b) } F_{mn } { } { \alpha\beta } 

6 //UTF8 format will print greek characters 

7 println t .toString (CrrFS) 

[> (a+b) ** (-1) *F_{mn} {a/3} 



CO . redbery . core . context . OuputFormat . * should be added to static imports. 
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8 //Wolf ramMathematica format 

9 println t .toString {Wolf ramMathematic4 

t> Power [b+a, -1 ] *F_{mn } { \alpha\beta } 



By default, Redberry uses its own output format, so line 3 is equivalent to 

println t . toString (J^edberry) 

This line produces string in Redberry input format, so the result can be parsed back into Redberry. 
The LaTeX format produces simple IM^i^X code that could be parsed by WI)^. If console 
supports Unicode characters, then the UTF8 format can be used to print Greek letters as is, 
while in other aspects this format is the same as default. Finally, the Wolf ramMathematica 
format produces the expression in form that can be directly parsed by Wolfram Mathematica. It 
is useful for further manipulations in Mathematica with symbolic (without indices) expressions. 

2.2 Tensors 

Redberry is written in Java and makes extensive use of its object-oriented features. Though 
user may not be familiar with object oriented programming (OOP) to use Redberry, it is still 
useful to understand primitive object types that are used in the CAS. There are three central 
object types in Redberry: Tensor, Indices and Transformation. Objects of the same class 
share common properties and can be manipulated in a common way^. Tensors are considered in 
this section, while indices and transformations are discussed in Sec. 2.5 and Sec. 3 respectively. 

Each mathematical expression in Redberry is a Tensor . Any Tensor have Indices and 
content (summands in case of sum, arguments in case of functions, etc.; thus, tensors in Redberry 
are containers of other tensors). Here is how these properties can be accessed via Groovy syntax: 



1 def product = 'F'^{A}_{B \\mu \\nu} * a'.t 

3 //getting element by index 

4 println product [ ] 

[> F''{A}_{B Wmu Wnu} 

5 //size of tensor as container 

6 //(here it means number of multipliers) 

7 println product . size ( ) 

[> 2 

8 //enumerating elements 

9 product . each { i -> println i } 

[> F''{A}_{B Wmu Wnu} 
> a 

10 //getting indices of some expression 

11 def sum = 'A'^ijk * B_i + N^j * M'^k'.t 

12 println sum. indices 

^In fact, Java internals of Redberry are highly introduced into main user interface (Groovy scripting). Despite 
the fact that only two developers created Redberry, the usage of Groovy as a general programming environment 
allowed us to implement rather convenient way for flexible Redberry usage. One of the main targets of Redberry 
is availability to the broader audience. So, authors tried to introduce simplicity and sameness as much as possible 
into the internal architecture as well as into Groovy binding. When architecting Redberry the authors were guided 
by intention to mimic the natural sense of tensorial expressions. 
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[> ''jk 



The meaning of the above hnes of code is pretty evident. First two examples demonstrate two 
ways of accessing (getting by index and enumerating) child elements of a tensor. The last example 
demonstrates that indices of a tensor can be taken via indices property. 

In order to create tensors programmatically, Redberry defines all arithmetic operations for 
Tensor objects (mathematical expressions)^: 



def t = 'A_i+C_i'.t, u = 'G'^i'.t, k = 'S_m'"m'.t 
println t*u+k 

[> (A_i+C_i) *G''i+S_m^m 

//one even can use this syntax to 
//create expressions from the scratch 
def X = 'x' .t, y = 'y' .t 

t = (-X + sin(y))**2 / (x + sin(-y)) + x - sin(y) 
println t 

[> 2*Sin [-y] +2*x 



However, such syntax may be inconvenient when dealing with indexed objects, since names of 
variables do not reflect structure of indices of expressions. 

Tensors are immutable and modification operations return the new instance: 



def t = 'A_i+C_i'.t 

def X = t.set(0, 'N_j*D'' j_i' .t) 

println x 

[> N_j*D^ j_i + C_i 

def y = t. remove (0) 
println y 

[> C_i 
println t 

[> A_i + C_i 



2.3 Einstein notation 

At this point the basic conventions arising from Einstein notation should be clarified. As 
was previously mentioned, Redberry distinguishes covariant (lower) and contravariant (upper) 
indices. The corresponding property of single index we call state. Two indices are considered 
to be contracted if and only if they have similar names and types but different states. As a 
consequence of this convention there are some natural restrictions on general structure of the 
expressions. 

First of all, the following notation (used, for example, in Euclidean space) is illegal in Red- 
berry: 

^The implementation uses Groovy operator overloading, while the Java equivalents are placed in Tensors 
class. 
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def t = 'F_aa' .t 

t> InconsistentlndicesException 



The correct input for the above tensor should be the fohowing: 
def t = 'F_a''a' .t 

In addition, the error occurs when expression is meaningless because of several indices with the 
same name and state coexist in the same product: 



//meaningless expression 
def t = 'F_i j*M'^i*N'' j*K'^ j' .t 

t> InconsistentlndicesException 



The entire architecture of Redberry built in such a way that the above illegal situations can not 
ever arise during manipulations and user do not need to take care about it. 

Other thing arises with dummy indices in products, where one or more multipliers are sums. 
It is convenient to write 

There is no problem here, since dummy index is in the scope of sum, which " outer indices are 
upper , while the first multiplier indices are lower ^j^. However, if we try to expand brackets 
in this expression naively (i.e. without relabeling of dummy index /i) we shall face the ambiguity 
as described in the previous examples. Thus it is better to relabel such dummy indices in sums 
right after parsing. Such transformation is automatically performed for every input expression 
in Redberry. Consider the following code: 



def t = 'F_mn* (A'^ab + M_m*N^mab) ' . t 
println t 

[> F_mn* (A^ab + M_c*N^cab) 



As we can see, dummy index m was automatically renamed to c , because conflict with free 
index of the whole tensor was detected. Similar behaviour is also valid for powers in products: 



def t = 'F_a* (A''a*B_a)**2' .t 
println t 

[> F_a* (A-^b*B_b) **2 



2.4 Standard form of mathematical expressions 

A core function of any CAS is its ability to reduce arbitrary expression to some standard form 
(SF), which is then used everywhere in manipulations. This approach facilitates comparison 
and matching of expressions and gives a way for more robust and fast algorithms of almost all 
transformations. Redberry uses the same paradigm, so any intermediate and resulting expression 
is guaranteed to be in the SF. 

Consider the following examples, which give an idea of SF in Redberry: 

^If consider sum as a one tensor T"" = A"!^ + M^N""'^ 
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println ' (a-b) +c+ (b-a) '. t 

[> c 

println ' c* (a-b) * (b-a) /c' . t 

[> - (a-b) **2 
println ' 2/3-27** (1/3) /9' . t 

[> 1/3 

println 'Sin[2 + 2 .*!]** (1/4) ' .t 

[> 1 . 38307 - . 144188*1 



The first line demonstrates that Redberry performs the reduction of similar terms. The second — 
that same (to within a sign) multipliers are collected into powers or reduced. Numbers are always 
collected and reduced if possible. If the expression contains floating-point numbers, then it will 
be completely reduced (calculated). This behaviour is similar to the majority of symbol-oriented 
computer algebra systems and needs no more detailed explanation. 

The additional conventions on standard form arises in expressions that contain tensors. The 
most remarkable convention is on the SF of sum. It is best to demonstrate it by example: 



println ' a*F_mn+ (a+b)*F_mn' . t 

[> (2*a+b) *F_{mn} 
println ' (x_a'^a+y_b^b)*X_m*X*m+ (z_n''n-y_d''d)*X_a*X''a' . t 

[> (x_{a}-^{a}+z_{n}-^{n}) *X-^ {m} *X_{m} 



As one can see, Redberry tries to factor out parts of products which contain all multipliers which 
have nonzero number of "outer" indices^. 

The other important detail is the ordering of summands and multipliers within sums and 
products. Elements of sums and products are sorted by their 32 bit hash codes. Hashes for 
simple tensors (e.g. x or k_p ) are generated randomly at each Redberry run, while the hashes 
of complex tensors (e.g. Sin [x] *k_i ) are calculated according to certain complex rules^. Usage 
of the pseudorandom generator allows to obtain nearly uniform distribution of hashes of tensors, 
which significantly improves performance. However, drawback of this approach is that ordering 
of expressions changes from run to run, and e.g. product a*b*c will be sorted differently 
at different runs (b*c*a or b*a*c etc.). Of course, all similar expressions will have similar 
ordering in the current session. Still, Redberry have tools to fix the seed of pseudorandom 
generator inside, so that expressions will have equal ordering from run to run. 

In contrast to many other tensors-oriented CASs, like xAct [5] or Cadabra [1], Redberry do not 
pay attention to the so-called indices canonicalization (sorting) problem. It uses fundamentally 
different approach for the problem of tensors comparison, which internally takes into account 
symmetry properties of tensors and does not depend on any particular canonical form of tensor 
indices. 

^Multipliers like (fi:''^ + u^a) are considered to have no "outer" indices, wile tensors like a;'"m are considered 
to have "outer" indices, because they contribute to a whole product indices. 

^The most important property of hash functions defined for complex tensors is its "insensibility" for particular 
names of indices but "sensibility" for their contractions. So, renaming of dummy or free indices is not affecting 
hash code, but changing the structure of contractions (e.g. contraction of two free indices) does. 
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In conclusion, authors want to emphasize that user can rely on fact that during any calculation 
any expression in Redberry is reduced to SF. This fairly simplifies implementation of custom 
transformations and algorithms. 

2.5 Indices 

Presence of indices of expressions is a main distinguishing property of tensor-oriented CASs. 
As mentioned above, there is an indices property defined for each expression in Redberry^. 
Returned object is an object of type Indices and have a number of methods and properties 
to work with it. Here are some examples of possible operations with indices objects^: 



1 def t = '2*x_am*f''m* (a'^n+b'Ti) '.t 

2 def ind = t. indices 

3 println ind 

[> ^{mn}_{am} 

4 println ind. size () 

[> 4 

5 println ind. free 

[> "{nl.ia} 

6 println ind. inverted 

[> ^{am}_{mn} 

7 println ind. upper as Indices 

> { mn } 



Methods used in the above example, are inherent in any indices object. Their names clearly 
implies their meaning'^. 

Although the presence of indices is inherent in all tensors, diff^erent types of expressions have 
different subtypes of Indices objects. This difference arises from availability or unavailability 
of information about indices order. Consider indices of simple tensor: 



def simple = ' F_{inn} { \\beta\\alpha}_{ba\\alpha} .'t 

println simple . indices 

[> _{mnba} ^ { \beta\alpha}_{ \alpha} 



We shall call indices of simple tensors as "simple indices". For "simple indices" the order of indices 
is defined, which is its main distinguishing property. In other words, permutation of indices will 
result in changing of mathematical sense of the particular simple tensor (unless this tensor is 
symmetric with respect to this permutation, see Sec. 2.6). However, as we previously mentioned, 

^ get Indices ( ) method in Java API 

^The last line in this example requires cc . redberry. core . indices . Indices to be added to im- 
ports. 

^FuU list of indices properties and methods as well as corresponding documentation can be found in the 
JavaDocs for the Indices class (see full documentation on Redberry web site). 
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each index belongs to some index type (e.g. Greek upper/lower case or Latin upper/lower case, 
etc.). Indices of different types are considered to have different mathematical nature (e.g. Greek 
indices are Lorentz, Latin are SU(N) etc.), so the relative position of indices with different 
types is not important. Thus, Redberry sorts indices of simple tensors according to their types, 
preserving the relative order of indices belonging to the same type (see the above example). 
Indices of TensorField are also '"simple indices" and follows the same ordering rules. 

Another type of indices is inherent in all other types of tensors. Consider the following 
product: 



def pr = 'F_{mn}*F'^{\\beta\\alpha)*F_{ba\\alpha} '.t 
println pr. indices 

[> { \alpha\beta }_{ abmnValpha } 



From the mathematical point of view, order of product indices is undefined. This allows to set 
up some certain ordering rules (mainly for technical reasons, related to performance). As one 
can see from the example, all indices are sorted according to the following rules: first upper then 
lower, first Greek then Latin, in each group of indices with the same type indices are sorted in 
lexical order. The similar rules are adopted for sums: 



def sum = ' R'^a_amn''\alpha+K^i_inm^\alphaJt 

println sum. indices 

[> { \alpha}_{nm} 



The only difference, is that according to the sense of sum, its indices are only free indices. 

All methods from the Listing 2.5 ( .inverted , .free etc.), return objects with the same 
rules of ordering as in the initial indices object. 

Single index At the low-level, Redberry stores each single index as 32-bit integer which en- 
codes all information about index: state, type and name (serial number in the alphabet). How- 
ever, there are many utility methods to operate with single index^, so, the user is free from 
knowing the concrete bits layout. In order to illustrate some possible manipulations with indices 
in Groovy scripting, consider the following advanced example^: 

Listing 2 : Get contractions between two tensors 

1 def getDummy = { 



2 tl, t2, indexType -> 

3 def indl = tl . indices . free, 

4 ind2 = t2 . indices . free, 

5 dummy = [ ] 

6 for (i in .. indl . size (indexType) - 1) 

7 for (j in .. ind2 . size (indexType) - 1) 

8 if (areContracted(indl [indexType, i], ind2 [ indexType, j])) 

9 dummy << getNameWithType (indl [indexType, i]) 
10 dummy as int [ ] 



11 } 

"""see IndicesUtils class 

^The following classes from cc . redberry . core . indices package should be added to imports: 
Indices , IndexType.* (statically) and IndicesUtils.* (statically). 
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12 def tl = 'T_{nm\\alpha}^{a\\beta}'.t 

13 def t2 = ' V'^{nc}_{a\\beta} '.t 

14 println getDummy(tl, t2, LatinLowsi^ as Indices 

[> _{an} 



In this example we define function which takes two tensors, index type and returns an array 
of contracted indices, having specified type, between specified tensors. The information about 
available index types is placed in the IndexType enum. Of course, the above code is not 
the most efficient implementation of such method, but aims to emphasize the ease of accessing 
low-level Redberry infrastructure. 



2.6 Symmetries of tensors 

The next distinctive feature of tensors is their symmetries. Consider symmetries under per- 
mutations of indices. Permutational symmetries in Redberry can be defined for indices of simple 
tensors and tensor fields. For example, in the following code, the symmetries are defined for Rie- 
mann tensor and then all possible permutations of indices (which follows from those specified) 
are enumerated: 



def t = 'R_abcd' .t 

addSymmetry {t , 2, 3, 0, 1) 
addAntiSymmetryit, 1, 0, 2, 3) 



//printing all possible permutations for this tensor 
for(s in t . indices . symmetries) 
println s 



[> 


[0, 


1, 


2, 


3] ( + ) 


[> 


[2, 


3, 


0, 


1] ( + ) 


[> 


[1, 


0, 


2, 


3] (-) 


[> 


[2, 


3, 


1, 


0] (-) 


[> 


[3, 


2, 


0, 


1] (-) 


[> 


[3, 


2, 


1, 


0] ( + ) 


[> 


[0, 


1, 


3, 


2] (-) 


[> 


[1, 


0, 


3, 


2] ( + ) 



//printing only basis 

symmetries . basis . each {s -> println s } 

[> [0, 1, 2, 3] ( + ) 
[> [2, 3, 0, 1] ( + ) 
> [1, 0, 2, 3] (-) 



As one can see, there are eight permutational symmetries (including identity), which can be 
obtained by combining the basic symmetries. Redberry supports permutational symmetries and 
anti-symmetries. Tensor with only one type of indices was used here for simplicity, but in case 
of several types of indices, symmetries can be set for each index type separately^. The example 
also shows method to access basis symmetries. 

Once set, symmetries of tensor affects all further manipulations with it. The following example 
is a good demonstration of this feature (code from the previous example assumed to be already 
executed) : 

^see JavaDocs for details 
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println ' R''abcc3*R_ef dc*R''ef_ab + R_rc''d&R_ab''rc*R_fci'^ba' . t 

[> 



Here zero was returned right after parsing. This is because Redberry automaticahy reduces sums 
to the standard form and 

R^'^'^'^RefdcR^^ ah = — Rrc^'^ Rab^'^ Rfd"" 

according to the specified symmetries. Such architecture requires user to set all symmetries of 
simple tensor before it will be parsed inside any complex structure like sum or product. If this 
criteria is not fulfilled, results of calculations may be erroneous. 

Another restriction concerns sequential addition of symmetries. Consider the following code: 



def t = 'R_abcd'.t 

addSymmetry (t. , 2, 3, 0, 1) 
addAntiSymmetry (t- , 1, 0, 2, 3) 
addAntiSymmetry (t- , 3, 2, 1, 0) 

[> InconsistentGeneratorsExcept ion 



Exception is thrown on the last line. It is easy to verify, that permutation used in this line is a 
combination of two previous, but sign is different. Thus, some combinations of symmetries and 
anti-symmetries can come into conflict, which causes an error in Redberry. 

Redberry provides tools to find permutational symmetries of complex tensors. Consider the 
following example^: 



addAntiSymmetry('R_abc', 1, 0, 2) 
addSymmetry {' K_a.h' , 1, 0) 

def t = ' (R_abc*A_de + R_bde*A_ac) ★A'^ce + R_adb'.t 

def symmetries = findIndicesSymmetries('_ah6.'.s±, t) 
for (s in symmetries) 
println s 

[> [0, 1, 2] ( + ) 
\> [2, 1, 0] (-) 



The first permutation is identity, while the second means that if Rabc = —Rbac and Aab = ^baj 
then tensor 

Tabd = {RabcAde + Rbde^acjA'^'^ + Radb 

satisfies the symmetry property Tahd = ~Tdba^- 

While the permutational symmetries of tensors are well covered in a number of tensors-oriented 
CASs, the so-called multi-terms symmetries like Bianchi identities, are faintly covered or absent 
at all in the majority of existing systems. In fact, authors know only one system — Cadabra [1], 
which fully supports multi-term symmetries. The basic idea utilized by Cadabra system is usage 

^ cc . redberry . core . utils . TensorUtils . * should be added to static imports. 

^Static method f indlndicesSymmetries (...) is defined in class TensorUtils . This method 
takes a simple indices as the first argument in order to specify the relative order of indices, for which one want to 
find symmetries. 
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of Young tableau projectors to reduce expressions to the simplified form. Consider the following 
identity^ : 

where Wabcd is a Weyl tensor. In order to proof this identity using a Young projector, one need 
to apply the following substitution to the above expression: 

Wabcd =\{2 Wabcd - Wadbc + W acbd) (2) 

This identity is derived from the Ricci cyclic identity and ordinary permutational symmetries of 
the Weyl tensor. The following Redberry code proofs the identity (1): 

Listing 3 : Multi-term symmetries 

1 addAntiSy/nmetry (' W_abcci', 0, 1, 3, 2) 

2 addSymmetry ('W_abcd', 2, 3, 0, 1) 

4 def t = { ' W''p_qT_s*W_p''t_r''u*W_tv^qw*W_u''vs_w' 

5 + ' - W^p_q''r_s*W_p^qtu*W_rvtw*W''sv_u'^w' 

6 + ' - W_mn'^ab*W'^n_pb''c*W'^ms_cd*W_s'^pd_a' 

7 + ' + l/4*W_mn'^ab*W'^ps_ba*W'^m_p^c_d*W'^n_s^d_c') .t 

9 def s = 'W_mnpq = 1 / 3* ( 2 *W_mnpq - W_inqnp + W_mpnq)'.t 

10 def r = (s & Expand) » t 

11 println r 

[> 



At the moment Redberry have no built-in functionality to construct substitutions like (2) based 
on the Young projectors. But, as could be seen from the above example, if one define the 
corresponding substitution manually, then Redberry allows to work with multi-term symmetries 
in the Cadabra way. Authors plan to incorporate support of Young projectors in the next release 
of Redberry. 

3 Transformations 

3.1 Applying and manipulating transformations 

All transformations in Redberry are first-class objects (which means that they can be assigned 
to variables) and share a common way for their applying and manipulating. Consider the syntax 
for applying transformation to mathematical expression: 



1 def tr = Expand 

2 def t = ' (A_k+B_k)*c' .t 

3 def rr = tr >> t, 

4 rl = t << tr 

5 assert rr == rl 

6 println rr 

^This example is taken from Section 2.2 of the 'Cadabra: reference guide and tutorial' by Kasper Peeters, 
which is available on Cadabra web site. 
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[> c*A_k+c*B_k 



As it seen from the example, transformations are applied using left shift << or right shift >> 
operator. It should be noted that both operators are left-associative, which means that, for 
example, A<<B<<C will be treated as "apply transformation B to tensor A, and then apply C to 
the resulting tensor", while C>>B>>A will be treated as "apply C to B and then apply resulting 
transformation to tensor A". The latter behaviour is not what one can expect, so it is better to 
use a special & operator, which allows to join a set of transformations into a single one: 



1 def t = ' (a+b)*c' .t 

2 def expandAndSubs = Expand & 'c = a + b'.t 

3 println expandAndSubs >> t 

[> a* (a+b) +b* (a+b) 

4 def subsAndExpand= 'c = a + b'.t & Expand 

5 println subsAndExpand >> t 

[> a* *2+2 *a*b+b* *2 



Transformations (like tensors) are immutable in Redberry. Some transformations may take 
required or optional arguments using square brackets: 



1 def elirtiinateWhileExpand = Expand[Elim±nateMetr±c4 r 

2 X = ' (g_mn + d_m'^a*g_na)*f'^mn' .t 

3 println elirtiinateWhileExpand >> x 

[> 2*f_m^m 

4 def diff = Differentiate! ' x_in' ] , 

5 y = ' x_m*x''in' . t 

6 println diff >> y 

[> 2*x^m 



In this example, the Expand transformation takes an optional parameter, which is a trans- 
formation to be applied on each level of expand procedure. In contrast, the argument of 
Differentiate transformation is required. In both cases a new object will be created and 
assigned to a corresponding variable. The meaning of arguments is specific for each particular 
transformation and will be discussed in further sections. 

3.2 Selected general-purpose transformations^ 
3.2.1 Substitutions 

The most frequent transformation in all computations is a substitution. Here we shall discuss 
the usage aspects of substitutions, while the idea of the underlying algorithms can be found in 
Sec. 5.1. 

^Transformations covered in this section are part of redberry-core , their Java implementation can be 
found in cc . redberry . core . transformations package. 
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The very important feature of any tensors-oriented CAS is automatic relabehing of dummy 
indices in the case of dummy indices clash. Redberry takes care about it in ah types of substitu- 
tions. Consider, for example, the following simple substitution: 

X = Xa in [xfa^ya){xh + Zh) 

Here is a code to perform this substitution in Redberry: 



def s = 'x = x_a''a'.t 

def t = ' (x*f_a + y_a)*(x*f_b + z_b)'.t 
println s >> t 

[> (x_{d}"{d}*f_{a}+y_{a}) * (x_{ c } " { c } * f_{b } +z_{b } ) 



As one can see, the appropriate relabelling was performed automatically. 

Redberry supports substitutions of tensor fields and automatically performs matching of the 
arguments of functions during substitution: 



def s = 'F_ij[x_m, y_m] = x_i*ry_j'.t 

def t = ' T^ab*F_ab [p^a - q'^a, p'^a + q'^a] '.t 

println s >> t 

[> T-{ab}* (p_{a}-q_{a}) * (p_{ b } +q_{b } ) 



If tensor field depends on indexed argument, then some ambiguities arises when mapping the 
arguments. For example, the following substitution: 

can be performed in two different ways: (a) matching Xij — )• XiUj gives XkVjf-', (b) matching 
Xij — )> XjUi gives UkXjf^. This is because the indices of product (xiyj) are not ordered by their 
nature. To explicitly specify the matching rule in such ambiguous situations, Redberry allows 
to enter the correspondence of indices in the field arguments: 



1 def s = 'F_i[x_mn] = x_i]»f^k'.t 

2 def t = 'F_k[x_i*y_j] ' .t 

3 println s >> t 

[> x_{ k } *y_{ a } *f '~ { a } 

4 t = ' F_k [x_i*y_j :_ji] ' . t 

5 println s >> t 

[> x_{ a } *y_{ k } *f ^ { a } 



The second line in this example is actually equivalent to 
def t = 'F_k[x_i*y_j:_ij] '.t 

because indices of product are sorted according to the rules discussed in Section 2.5. 

As well as simple tensor and tensor field substitutions, Redberry fully supports all other types 
of substitutions, taking into account both indices symmetries and indices contractions. Consider 
the following complex example. Let's apply the following substitution 
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to tensor 

where Rabc is antisymmetric: Rabc = ~Rcba- It is clear that the sum over first, third and fom'th 
terms can be rewritten in the following way: 

fi - Rhj + RijkF^^ = fi + Rj,k F''' - RkjiF^' = n + Rua F'^ - RijiF^' , 

where the last step is a simple dummy relabelling. After this steps, it becomes obvious that these 
terms match the l.h.s. of the substitution. In Redberry one can e.g. do 



addAntiSymmetry {' R_mnp' , 2, 1, 0) 

def s = ' f _m + R_bma*F''ba - R_ljmtF^lj = R_bai»F''ab' . t 

def t = ' f _i + R_ijk*F''jk + R_ijl»F''kj - R_ki j*F'' jk' . t 
println s >> t 

[> 



We see, that Redberry matched the l.h.s. of the substitution in tensor t and automatically 
reduced the resulting sum to the standard form, which, in turn, gives zero. 

Redberry takes into account not only predefined symmetries of simple tensors, but also sym- 
metries of any complex expression, which arise from its structure. For example: 



def s = 'K_a * (A'^ab - A^ba) = F''atA_a'^b' . t 

def t = 'K_p * (A'^qp - A'^pq) + F''btA_b^q' . t 

println s >> t 

[> 



The result is zero since tensor {Act — Af,c) is antisymmetric. As well, Redberry takes care about 
symmetries of indexless objects: 



def c = 'Cos[a - b] = c'.t, 
s = 'Sin[a - b] = s'.t, 

t = 'x = Cos[b - a]**3 + Sin[b - a]**3'.t 
println (c | s) >> t 

[> X = c**3 - s**3 



As one can see, the built-in definitions of sin and cos are set up to be odd and even respectively. 

There is an important note on applying several substitutions at a time using the joining of 
transformations. Substitutions joined with & operator will be applied sequentially. However, 
sometimes it is necessary to apply several substitution rules "simultaneously". Consider the 
following example: 



def X2YandY2X = 'x=y'.t & 'y=x'.t, 
X2YorY2X = 'x=y'.t | 'y=x'.t, 
tensor = 'x+2*y'.t 

println X2YandY2X >> tensor 
t> 3*x 

println X2YorY2X >> tensor 
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[> y+2*x 



The first transformation ( x2YandY2x ) means just sequential applying of two provided substitu- 
tions, while the second ( x2YorY2x ) performs both substitutions "simultaneously". 
The last thing need to be discussed in this section is that the expression like 

def t='x=a+b'.t 

is a mathematical expression ( Tensor ) and transformation ( Transformation ) at the same 
time. In order to construct such expression programmatically from a given l.h.s. and r.h.s. one 
can use the following syntax: 



def Ihs = 'x'.t, rhs = 'a + b'.t, 

subs = Ihs.eq rhs 
println subs >> 'x**2'.t 

[> (a + b) **2 



3.2.2 DifTerentiate 

This transformation allows to take derivatives with respect to indexed objects. Consider the 
examples: 



1 def tensor = ' Sin [f_ab*f '^ab] ' . t 

2 println Dif ferentiate[ ' f_mn' ] >> tensor 

[> 2*Cos[f"{ab}*f_{ab}]*f"{mn} 

3 //derivative of antisymmetric tensor 

4 setAntiSymmetric {' R_ah' ) 

5 println Differentiate[ ' R_ab' ] >> 'R_mn'.t 

[> (1/2) * (d_{m}"{a}*d_{n}'-{b} - d_{ n } " { a } *d_{m} " {b } ) 



The Differentiate transformation takes care about dummies relabelling and symmetries 
of tensors. The following convention is adopted: 



where is a number of elements in the sum, and r.h.s of this expression have the same symmetries 
as l.h.s. 

It is also possible to pass additional transformations which should be applied after each step 
of differentiation (mainly for performance reasons): 



1 //setting up symmetries of Riemann tensor 

2 addAntiSymmetry('R_abcd', 1, 0, 2, 3) 

3 addSymmetryC R_ahcd' , 2, 3, 0, 1) 

5 def tensor = ' R''acbd*Sin [R_abcd*:R''abcd] ' . t, 

6 varl = ' R''ina_m'^b ' . t , 
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var2 = 'R^mc_m^d'.t 



9 def diffl, diff2 

10 timing { 

11 //take second derivative and 

12 diffl = {DiffBrentiateivar2, 

tensor } 

[> Time : 1338 ms . 

13 timing { 

14 //take second derivative and simplify permanently 

15 diff2 = Diffsrentiatdivarl, varl, ExpandAndEliminat^ >> tensor } 
[> Time: 14 ms . 

16 assert diffl == diff2 



As one can see from this example, if ExpandAndEliminate (which simply expand out brackets 
and eliminates contractions with metrics and Kronecker deltas) performs on the each level of the 
differentiation, the derivative will be taken much faster, then if it applies only to the result. 

3.2.3 EliminateMetrics 

This transformation eliminates metric tensors and Kronecker deltas, which are contracted 
with other tensors: 



1 def tensor = ' g_nm*A''m*d''n_a' . t 

2 println EliminatsMetrics » tensor 

> A_a 

3 tensor = ( ' g''mn*g''ab*g''gd* (p_g*g_ba + p_a*g_bg) ' + 

4 ' * (p_m*g_dn + p_n*g_dm) ' ) . t 

5 //eliminate metrics in D dimensions 

6 println { {EliminateMetrics & 'd_a'^a = D'.t) >> tensor) 

[> 2* (1+D) *p'^ {d} *p_{d} 



3.2.4 Expand 

There are several transformations which expand out products and integer powers of sums: 
Expand , ExpandAll , ExpandNumerator and ExpandDenominator . These well-known 
transformations needs no special explanation except the fact that they can take an additional 
transformations as an arguments to be applied on the each level of expand procedure. Consider 
the following example: 



1 def tensor = ( ' (g_af*g_bc+g_bf*g_ac+g_cf*g_ba) * (T_d*T_e+g_de) ' + 

2 ' * (g''db*g''ae + g'^de*g'^ab) ' ) . t 

4 //eliminates metrics in four dimensions 

5 def eliminate = EliminateMetrics & 'd_a^a = 4'.t 



then simplify 

varl] & ExpandAndEllmlnat4 » 
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7 //expand and then eliminate 

8 def rl = {Expand & eliminate) >> tensor 

9 println rl 

[> 2*T_{c}*T_{f }+30*g_{fc}+7*g_{fc}*T_{e}*T"{e} 

10 //eliminate while expand 

11 def r2 = Expand[ eliminate] >> tensor 

12 assert rl == r2 



As we see, lines 6 and 8 produce the same result, but the latter spends less time in the calculation 
because it applies additional simplifications not only to the final result, but to all intermediate 
tensors, drastically reducing their complexity. 

As was many times noted, all operations in Redberry are safe from conflicts of dummies, and, 
for example, Expand transformation automatically renames dummy indices of integer powers 
in the resulting sum: 



1 //dummy indices in power 

2 def t = ' (a+b) / (A_m^m+B)**2*A_a' .t 

3 println ExpandDenominator » t 

[> (a+b) * (2*B*A_{b} {b} +A_{m} ^ {m} *A_{b} {b}+B**2)**(-l) *A_{a} 



3.2.5 Factor 

This transformation factors symbolic polynomials over the integers in expressions^. Consider 
the examples: 



1 //univariate polynomial 

2 def t = '1 + 2*x + x**2' .t 

3 println Factor >> t 

> (l+x)**2 

4 //multivariate polynomial 

5 t = '2*x**3*y - 2*a**2*x*y - 3*a**2*x**2 + 3*a**4'.t 

6 println Factor » t 

t> (x+a) * (x-a) * (-3*a**2 + 2*y*x) 

7 //tensorial expression with symbolic polynomial parts 

8 t = Expand » ' (a+b)**4*F_mn + (x**6-y**6) *R_mn' . t 

9 println Factor » t 

t> (a+b) ** 4 *F_mn+ (x+y ) * (x-y ) * (x*y+x* *2+y**2 ) * (-x*y+x**2+y**2) *R_mn 



Factor traverses the expression from head to children and factors all symbolic (indexless) 
polynomials (both univariate and multivariate) and rational functions. So, in the case of com- 
pletely symbolic expression it applies only to the top algebraic level. In the case of rational 

^The code of polynomials factorization is taken from free and open source Java Algebra System [6, 7] designed 
by Heinz Kredel. 
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expression, Factor first calls Together (Sec. 3.2.6), then factors numerator and denomina- 
tor. 

It should be noted, that in the current version of Redberry, performance of Factor transfor- 
mation in case of multivariate polynomials may be low in some cases and some huge expressions 
may not be factored completely. Also Factor applies only to expressions, which does not con- 
tain trigonometric functions, tensor fields or non integer powers of variables. These issues will 
be fixed in the upcoming release of Redberry. 

3.2.6 Together 

There are two transformations, which put terms in a sum over a common denominator: 
Together and TogetherFactor . The last one also cancels symbolic (without any indices) 
factors in the result. Consider the examples: 



1 def t = 'a/b + c/d'.t 

2 println Together » t 

> b** (-1) *d** (-1) * (b*c+a*d) 

3 //to common denominator and do not cancel factors 

4 t = 'x**2/(x**2 - 1) + x/(x**2 - l)'.t 

5 println Together >> t 

t> (-l+x**2) ** (-1) * (x+x**2) 

6 //to common denominator and cancel factors 

7 println TogetherFactor >> t 

> (1+x) ** (-1) * (-1+x) ** (-1) * (x+x**2) 

8 //together with tensorial expressions 

9 t = 'f_m/a + k_in/ (f.Mrfm) ' .t 
10 println Together » t 

[> a**(-l)*(f_{a}*f"{a})**(-l)* { f_{b } *f " { b } *f_{m} +a*k_{m} ) 



As it seen from the last example. Together effectively relabel dummy indices when it is neces- 
sary. 

3.3 Selected physical transformations^ 
3.3.1 DiracTrace 

This transformation calculates the Dirac trace of expression in four dimensions. First, consider 
traces without 75^: 



1 //set up matrices 

2 defineMatrix 'G_a', Matrixl. matrix 

3 //dirac trace transformation 

4 def dTrace = Di racTr a ce[ ' G_a' ] 



Transformations covered in this section are part of redberry-physics . The Java implementation can 
be found in cc . redberry . physics . feyncalc package. 

^ cc . redberry . core . indices . IndexType . * should be added to static imports 
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6 //trace of product of matrices 

7 def t = 'Tr [G_a*G_b] ' .t 

8 println t 

[> G.a'^a'.b' *G_b''b'_a' 

9 //calculate trace 

10 println dTrace >> t 

[> 4*g_ab 

11 t = 'Tr[(p_a*G^a + m)*G_m* (q_a*G''a-m)*G_n] ' .t 

12 println dTrace >> t 

[> 4 *p_{m} *q_{ n } +4 *p_{ n } *q_{m} -4 *m* *2 *g_{mn } -4 *p^ { a } *g_{mn } *q_{ a } 

13 t = 'Tr[G_a*G_b*G_m*G''a*G_n*G_f*G''b*G''f] ' .t 

14 println dTrace >> t 

[> - 3 2 * g_mn 



The first line tells Redberry to consider tensor G_a as matrix with additional one upper and one 
lower index of type Matrixl (Latin lower letters with strokes). This means that everywhere 
at the input Redberry will consider G_a as matrix G_a^i'_j' and use matrix multiplication 
rules for the products of such tensors (see Appx. A for details). Line 4 defines DiracTrace 
transformation with specified notation for gamma matrix. 
Consider traces with 75: 



1 //set up matrices 

2 defineMatrices 'G_a', 'G5', Afa trixl. matrix 

3 //setting up symmetries of Levi-Civita 

4 setAntiSymmetric {' e_abcd') 

5 //dirac trace transformation 

6 def dTrace = Di ra c Tra ce[ ' G_a' , 'G5', 'e_abcd'] 

8 def t = 'Tr [G_a*G_b*G_c*G_d*G5] ' .t 

9 println dTrace >> t 

O -4*I*e_{abcd} 

10 t = 'Tr[G_a*G_b*G_c*G''d*G_e*G_f*G5] ' .t 

11 println dTrace >> t 

[> 4 *I*g_ab*e_f c''d_e + 4*I*d_c^d*e_f abe+4 *I*g_fb*e_c^d_ea-4 *I*g_ce* 
e_f ab^d+4*I*d^ d_e * e_f abc-4*I* g_f a * e_c '~ d_eb 

12 t = 'Tr[(p_a*G^a + m)*G_m*G5* (q_a*G''a-m)*G_n] ' .t 

13 println dTrace >> t 

[> -4*I*p_{b} *q_{a} *e^ {a}_{n} {b}_{m} 



As we see from line 6, when expression contains 75, it is also necessary to specify the notation 
of 75 ( G5 ) and Levi-Civita tensor ( e_abcd ) to DiracTrace . 

3.3.2 UnitaryTrace 

This transformation calculates traces of unitary matrices: 
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1 //set up matrices 

2 defineMatrix ' T_A' , Mat rix2. matrix 

4 //structure constants are antisymmetric 

5 setAntiSymmetric (' f_ABC' ) 

6 //d-constants are symmetric 

7 setSymmetric {' d_ABC') 

9 //unitary trace transformation 
10 def uTrace = UnltaryTracdi ' T_A' , 'f_ABC', 'd_ABC', 'N'] 

12 //trace of product of matrices 

13 def t = 'Tr [T_A*T_B] ' .t 

14 println t 

[> T_A^A'_B' *T_B'^B'_A' 

15 //calculate trace 

16 println uTrace >> t 

[> (l/2)*g_AB 

17 t = 'Tr[T_A*T_B*T_C] ' .t 

18 println uTrace >> t 

[> (1/4*1) *f_{CAB}+ (1/4) *d_{CAB} 

19 t = 'Tr[T_A*T'^A + l]'.t 

20 println uTrace >> t 

[> N-1/2+ (1/2) *N**2 

21 t = 'Tr[T_A*T_B*T_C*T'^A] ' .t 

22 println uTrace >> t 

[> ( - (1/4) *N** (-1) + (1/4)*N) *g_{BC} 



As one can see, it is necessary to specify the notation used for structure and d- constants ( f_ABCD 
and d_ABCD ) and dimension of unitary group ( N ). 

3.3.3 LeviCivitaSimplify 

This transformation simpHfies combinations of Levi-Civita tensors: 



1 //three dimensions 

2 setAntiSymmetric {' e_abc' ) 

3 def t = 'e_abc*e^abd' .t 

5 //simplify in Euclidean space in three dimensions 

6 println IieviCivitaSimpIify euclidean [' e_abc' ] >> t 

[> 2*d"d_c 

7 //simplify in Minkowski space in three dimensions 

8 println IieviCivitaSimpIifyminkowski [' e_abc' ] >> t 

[> 2*d"d_c 
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10 
11 



t = ' e_abc*e^abc' . t 

//simplify in Minkowski space in three dimensions 
println IieviCivitaSimpIifyminkowski [ ' e_abc' ] >> t 



> 6 



12 //four dimensions 

13 setAntiSymmetrici' e_abcd') 

14 t = ' 4*I*e''h_d'^fb*e_abch*e_e^d_gf '.t 

15 //simplify in Euclidean space in four dimensions 

16 println IieviCivitaSimpIify euclidean [ ' e_abcd' ] >> t 

> 16*I*e_{eagc} 

17 t = ' 4*I*e''h_d'^fb*e_abch*e_e^d_gf ' .t 

18 //simplify in Minkowski space in four dimensions 

19 println IieviCivitaSimpIifyminkowski [ ' e_abcd' ] >> t 

[> -16*I*e_{eagc} 



As one can see, it is necessary to specify whether space is considered to be Euchdean or Minkowski 
and the notation for Levi-Civita tensor. The difference between Euchdean and Minkowski Levi- 
Civita tensors appears for even number of space-time dimensions (compare hues 16 and 19). Also, 
as it is seen for example from line 11, the transformation automatically substitute the dimension 
of space, which is considered to be equal to the number of Levi-Civita indices. 

4 Selected physical applications 

4.1 Feynman diagrams 

One of the most prominent applications of Redberry is calculation of Feynman diagrams 
in quantum field theory. As previously mentioned, Redberry provides several common physical 
transformations such as traces of Dirac gammas and SU(N) matrices, simplification of Levi-Civita 
combinations etc. It also provides a powerful tools for simple inputting of noncommutative matrix 
expressions. In this section as an example we shall consider the calculation of the well-known 
differential cross section of the Compton scattering in quantum electrodynamics (i.e. spinor 
electrodynamics, in contrast to the scalar electrodynamics, which was illustrated in Sec. 1.2). 
The physical background can be found in Sec. 5.5 of [8]. 

First of all it is necessary to note, that Redberry does not support noncommutative products of 
indexless variables. However, many noncommutative objects, which occur in physical calculations 
(e.g. spinors or Gamma matrices) have a matrix origin, which means, that they have additional 
matrix indices, which are usually omitted for convenience. For example, when we write product 
of Dirac bispinors and gammas, like e.g. 



we mean, that the above quantities have additional matrix indices, which can be written explic- 
itly: 



In this expression a',b',c' are special matrix indices. From such point of view these "noncom- 
mutative" products can be represented as ordinary products of indexed objects. Such matrix 
indices should have a nonmetric type, which implies that it is impossible to perform raising or 
lowering or define a symmetry which mixes indices with different states. Redberry provides an 
internal facilities allowing to input matrix expressions in a convenient form like (3). The matrix 




(3) 




(4) 
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indices will be automatically inserted at parsing and subsequent processing could be made as 
with usual indexed expressions. Detailed description of this feature can be found in Appendix A. 

So, let us turn to the physical aspects of the problem. The following lines give squared matrix 
element of the Compton scattering in quantum electrodynamics^: 

Listing 4 : Compton scattering in QED 

1 defineMatrices 'G_a', 'V_i', 'D[x_m]', Matrixl. matrix, 

2 'vu[p_a]', Matrixl. vector, 

3 'cu[p_a]', Matrixl. covector 

5 //photon-electron vertex 

6 def V = 'V_m = -I*e*G_m' .t, 

7 //electron propagator 

8 D = 'D[p_in] = -I* (m + p_m*G^m) / (m**2 - p_in*p''m) ' . t , 

9 //diagram a) 

10 Ma= ' cu [p2_m]*V_m*e''m[k2_m]*D [kl_m+pl_m]*V_n*e^n [kl_m]*vu [pl_m] '.t, 

11 / /diagram b) 

12 Mb= ' cu [p2_m]*V_m*e''m[kl_m]*D [pl_m-k2_m]*V_n*e^n [k2_m]*vu [pl_m] '.t, 

13 //matrix element 

14 M = Ma + Mb 

15 M = (V & D) >> M 

17 def mandelstam = setMandelstam{ 

18 ['pl_m': 'm', 'kl_m': '0', 'p2_in': 'm', ' k2_in' : '0']) 

19 M = (ExpandAll & mandelstam) >> M 

21 //complex conjugation 

22 def MC = M 

23 MC = 'vu[pl_m]*cu[p2_m] = vu [p2_m>cu [pl_m] ' . t >> MC 

24 MC = (Conjugate & InverseOrderOfMatrice^Matrixl] ) » MC 

26 //squared matrix element 

27 def M2 = ExpandAll » (M * MC / 4) 

29 //photon polarizations 

30 M2 = 'e_m[kl_a]*e_n[kl_a] = -g_mn'.t >> M2 

31 M2 = 'e_m[k2_a]*e_n[k2_a] = -g_mn'.t » M2 

33 //electron polarizations 

34 M2 = ' vu[p2_m]*cu[p2_m] = m + p2''nirG_m' . t >> M2 

35 M2 = 'vu[pl_m]*cu[pl_m] = m + pl^rarG_m' . t >> M2 

37 //applying trace of gamma matrices 

38 M2 = DiracTracei ' G_a' ] » M2 

40 //final simplifications 

41 M2 = {ExpandAndEliminate & ' d''m_m = 4'.t & mandelstam) >> M2 

42 M2 = 'u = 2*m**2 - s - t'.t » M2 

43 M2 = Factor » M2 

45 println M2 



cc . redberry . core . indices . IndexType . * should be added to static imports 
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[> 2* (-m**2+s+t) ** (-2) *e**4* (-8*s**2*m**2*t+4*s**3*t+2*s**4+t**3*s 
+2*m**8 + 4*m**4* s *t-m* *2 *t * * 3-2 *m* *2 *t * *2 *s + 3 *m* *4*t**2-8*s 
* * 3 *m* *2 + 12 *s * *2 *m* * 4 + 3 *s * *2 *t * *2-8 *m* * 6*s ) * (m**2-s) ** (-2) 



The above code reproduces the standard steps of Feynman diagrams calculation and prints 
squared matrix element of the Compton scattering averaged over polarizations of initial particles 
and summed over polarizations of final particles: 

(m^ — s)^(— + s + t)^ 

X (-Ss^m^i + is^t + 2s'' + t^s + 2m^ + im^st - m^t^ 

-2m^t^s + 3mH^ - 8s^m^ + Us'^m'^ + 3sH^ - 8m^s) 

Let's consider the code in more detail. Firstly, in lines 1 — 3 we tell Redberry to consider 
some input objects as (noncommutative) matrices, vectors or covectors. These objects are: 
gamma matrices ( G_a ), photon-electron vertex ( v_i ), electron propagator ( D [p_m] ), elec- 
tron wave function u{p) ( vu [p_m] ) and its conjugation u{p) ( cu [p_m] ). To make this, 

def ineMatrices (...) function is executed with lists of tensors followed by target object 
descriptors. Each descriptor specifies the type of matrix indices and a type of matrix (e.g. 

Matrixl .matrix , Matrixl . covector ). Here only one matrix index type is used. As result, 
for example, the combination 7^ ti(p2) will be treated as scalar with respect to matrix 

indices: 



def t = ' cu [pl_m]*G_m*vu [p2_m] ' . t 
println t 

[> cu_{a' } [pl_rti] *G_{m}"{a' }_{b' }*vu"{b' } [p2_m] 
println t . indices . free 

[> _{m} 



while the combination u{pi)u[p2) will be treated as matrix: 



def t = 'vu[pl_m]*cu[p2_m] '.t 
println t 

\> vu_{a' } [pl_m] *cu^ {b' } [p2_m] 
println t . indices . free 

[> "{b' }_{a' } 



The next step is to build the matrix element. There are two Feynman diagrams corresponding 
to the Compton scattering in the leading order: 



spins 
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This done in the hnes 5 - 14, where e_n[k_m] denotes photon polarization. Important point 
can be seen in the hne number 8, where the propagator is introduced. The identity matrix 
(Kronecker delta) is inserted automatically in the numerator, where the gamma matrix is added 
to the scalar mass: 



println 'm + p_i*G''i'.t 

[> m*d"{a' }_{b' } + p_i*G"{i a' }_{b' } 



At the next step (lines 17 - 19) we define mandelstam variables and expand out the matrix 
element. Let us consider the function setMandelstam { . . . ) : 

setMandelstami [' pi' :' ml' , 'p2':'m2', 'p3':'m3', 'p4':'m4']) 

It takes a map 'momentum — mass of particle' as an argument and returns a collection of 
substitutions, followed from the rules: 



So, for example, the line 17 will produce the following substitutions: 



println mandelstam 

t> pl_{m} *pl'^ {m} = m**2 

[> kl_{m} *kl'^ {m} = 

[> p2_{m} *p2'' {m} = m**2 

[> k2_{m} *k2'^ {m} = 

[> 2*kl'^ {m} *pl_{m} = s-m**2 

[> 2*k2'^ {m} *p2_{m} = s-m**2 

[> -2*p2^ {m} *pl_{m} = -2*m**2+t 

[> -2*kl_{m} *k2'^ {m} = t 

t> -2*k2'^ {m} *pl_{m} = u-m**2 

[> -2*kl_{m} *p2'" {m} = u-m**2 



which then will be applied in further calculations. 

The next step (lines 22 - 24) is to make up a complex conjugation of the matrix element. As 
it is known from physics, this can be done using the following relation: 



u 



s 



t 



{Pi - P-if 
{Pi - Pif 



{P3 + Pif 
{P2 - Pif 
{P2 - P3f 
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where 

M = ^ Ci 0.10.2 ... Oi 
i 

M = c* Oi Oi-i . . . oi 

i 

and usual notation a = J^a^ is used. This transformation is equivalent to two substitutions: 

u{pi) u{p2) u{p2) u{pi) 
M ^ M 

The first substitution is applied in the line number 23 and the second in the line number 24. The 
meaning of inverseOrderOfMatrices [Matrixl ] transformation is clear from its name: 



defineMatrix 'A_\\inu', Matrixl. matrix 
def p = ' A_\\mu*A_\\nu' .t 
println p 

[> A.iWmul'ia' }_{b' } *A_{ \\nu} " {b' }_{c' } 
println InverseOrderOfMatrice^Matrixl] >> p 

[> A_{\\nu}"{a' }_{b' l^A.iWmul'fb' }_{c' } 



At this point we can define squared matrix element in the line number 27 and perform the 
summation over photons (lines 30 - 31) and electrons (lines 34 - 35) polarizations. On the last 
step all products of gamma matrices are automatically converted to combinations of their traces. 
So the next step is to apply DiracTrace transformation. After this, it is necessary to eliminate 
metric tensors and apply mandelstam substitutions (lines 41 - 42). This step will produce a 
symbolic expression, which then can be finally simplified using the Factor transformation. 
The whole calculation takes less than 0.5 seconds on the hot Java Virtual Machine (on Intel 
Core 15 @ 2.27 Ghz). 

4.2 One- loop counterterms of arbitrary Lagrangians 

Another remarkable application of Redberry is calculation of the divergent part of the one- 
loop effective action of arbitrary Lagrangians. The theoretical formalism based on the extended 
t'Hooft and Veltman method of the background calculations was developed in [9] and successfully 
applied for a number of theories using REDUCE computer algebra system [10]. This algorithm 
calculates the one-loop counterterms for an arbitrary theory and background in four dimensions 
in curved space-time in the dimensional and higher derivative [11] regularizations. Redberry 
implements the algorithms developed in [9] for second and fourth order differential operators 
and provides a simple and convenient user interface^ . 

Let's start with a minimal physical background. It is well-known, that one-loop effective 
action for a general field theory with a given action S[(j)] can be expressed in terms of derivative 
of action with respect to the fields: 

r«=fnTrfln-^V 
^The implementation is placed in cc . redberry . phys ics . oneloopdiv package. 
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where (fii denotes a "background field" and Latin letters denotes the whole set of its indices"^. So, 
the main quantity, which determine the effective action is a differential operator 

In the most general case, this operator has the following form: 







I V V ^2 • • ■ 




+ 


glJ.llJ.2-- 






+ 




■-^lL-2 J y v 

I V V ^2 


■■■Vm^_2 


+ 


]\J-tJltJ2- 


t V jitj V jU2 • 




+ 




■■■^|.L-4 J v v 


• • • 



+..., (6) 

where is a covariant derivative with respect to space-time and gauge indices: 

where T^-y is a Christoffcl symbol and uj^i^ is a connection on the principal bundle. Commuting 
covariant derivatives it is always possible to make tensors K, S, W, N, M symmetric in the Greek 
indexes and we shall assume this condition in the further reading. 
It also required to introduce the following quantities: 

{Kn)] = K^^f^^-t^\^ n^^n^, ...n^^, (7) 
{Kn)-\^Kn)/ = 5,^ (8) 

where is a unit vector. The second equation defines tensor {Kn)~^ inverse to tensor Kn, 
which is an input tensor for the algorithm^. The other required input is a curvature tensor with 
respect to the principal bundle: 

[V^,VA^i = F^J^j (9) 

Given a set of tensors K, W , M, {Kn)~^ and F as input data, Redberry allows to calculate 
counterterms for the general second and fourth order operators of the following form 

L>(2)i = i^^VVMV. + W (10) 
L>(4) i = V^V.VaV^ + W^i^ V^V. + (11) 

Let us illustrate the usage by the particular examples. 



4.2.1 Non-minimal vector field 

Quantum theory of non-minimal vector field in a curved space-time is determined from 
the following action: 

^Hcrc the actual type of these indices is not important (e.g. this set can contain both space-time and SU(N) 
indices), so Latin letters used only for convenience. 

^Redberry also provides facilities to find out the inverse tensor from the given equation. See 
cc . redjaerry-. physics . utils . I»,ver?eTe»SQr class for examples. 
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where tensor P^^^ is symmetric. One-loop counterterms for this theory were found in [12] and 
[13]. In the case of vector field, the field indices marked with Latin letters in (5) are space-time 
indices, so the curvature tensor introduced in (9) is equivalent to Riemann tensor: 

The differential operator (5) corresponding to the vector field is: 

where □ = g^'^'V and A = 1 + 1/^. This is a second order operator, and in order to rewrite it 
in the form (10), it is necessary to symmetrize the second term by commutation of the covariant 
derivatives: 



/3 
a ) 



where Raj3 is the Ricci tensor. 

Finally, using equations (7) and (8) it can be easily found that 

{Kn)-\P = 5i + Y^n^n^ 

Hereby, at this point we have the whole set of the inputting data which is required by the 
algorithm: 

F ^ - B ^ 

K>^-J = g^^-5i - ^ [g^^^S':, + g'^^S^) 
Sf'J = 

t^/ = + ^RJ 
{Kn)-\^ = 5f + Y^n^n^ 

In further calculations we shall use the definition A = 7/(1-1-7) for convenience (so 7 = A/(l — A)). 

The following code calculates one-loop counterterms of the vector field theory in curved space- 
time (here g used for 7): 

Listing 5 : One-loop counterterms of the theory of vector field. 

1 setSym/netric ( ' P_\\inu\\nu') 

3 def KINV = ('KINV_\\alpha'^\\beta = ' + 

4 'd_\\alpha^\\beta + gni_\\alpha*n'^\\beta') .t 

6 def K = {'K^\\mu\\nu_\\alpha''\\beta = ' + 

7 'g''\\mu\\ni>d_\\alpha'^\\beta - g/ (a (1+g) )* (' + 

8 'g'^\\mu\\beta*d_\\alpha^\\nu + ' + 

9 'g''\\nu\\beta*d_\\alpha^\\mu) ') .t 

11 def S = 'S''\\rho'^\\mu_\\nu = O'.t 
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13 def W = {'W^\\alpha_\\beta = ' + 

14 'P^\\alpha_\\beta+g/ (a (1+g) )*R'^\\alpha_\\beta') .t 



16 def F = 'F_\\mu\\nu\\alpha\\beta = R_\\mu\\nu\\alpha\\betat 

18 //calculates one-loop counterterms of the second order operator 

19 def div = oneloopdiv2 {KINV , K, S, W, F) 



21 def counterterms = Factor >> div . counterterms 

22 counterterms = ' P^\\alpha_\\alpha = P'.t >> counterterms 

23 println counterterms 



t> counterterms = (1/120) * (-32 + 5*g**2 + 10*g) *R'' { \epsilon}_{ \mu} *R_ 
{ \epsilon} ^ { \mu}+ (1/48) *g**2*P**2+ (1/240) *R**2* (2 8+5*g 
**2+2 0*g) + (l/24)*(g**2 + 12 + 6*g)*P''{ \beta}_{ \alpha} { \alpha} 
_{\beta}+ (1/12) *g* (4+g) *R_{ \nu } { \epsilon } ^tP-^ { \nu }_{ \epsilon 
} + (1/24) *R* (g**2+4+2*g) *P 



In order to obtain one-loop counterterms in the dimensional regularization, one should multiply 
the result produced by Redberry by l/167r((i — 4) and integrate it over the space-time volume: 

= [ d^x^ f — (-32 + + 107)i?eui?''' + — 7^P^+ 

167r(d-4)7 ^ ^V120 ^ 48' 

+ ^^'(28 + + 2O7) + ^(7' + 12 + 67)P/3aP"^+ 

The above code is clear enough, but some remarks are needed. First of all, the main method 
oneloopdiv2 (...) , which calculates the counterterms of the second order operator, returns a 
special object, which holds some intermediate results (like e.g. RR, RF, FF parts of effective 
action from [9]). The whole result can be obtained by getting the value of .counterterms 
property. All input expressions must be in the same notation as in the original work [9] except 
tensor {Kn)~^, which should be denoted as KINV . At this moment, the implementation requires 
that all indices of input tensors should be lower Greek indices. Also, it assumed that field indices 
are placed at the end, so, for example, the first two indices of tensor K^^"^ p are contracted with 
covariant derivatives in (6), while the last two indices corresponds to the indices of vector field. 
Redberry does not support nonzero tensor S from (6), however it should be specified explicitly, 
like it is done in line 11. 

4.2.2 Minimal fourth order operator 

As an example of the fourth order calculation, let us consider minimal fourth order operator 
of the form 

Di^ = Vn' + W^^'iN^V, + M,^ (12) 
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After symmetrizatioii of the first term, we obtain the following input: 



K^^'^'J = (1/3) 5 J ig^^g'^' + g'^^g^' + g-^'g^ 



NPj = 

Ma" = 




The following code calculates one-loop counterterms for this operator: 

Listing 6 : One-loop counterterms for the minimal fourthe order operator 

1 def KINV = 'KINV_\\alpha''\\beta = d.Walpha-^Wbeta.'t 

2 def K = ('K'^\\mu\\nu\\gamma\\delta_\\alpha^\\beta = ' 



8 def S = 'S'^\\mu\\nu\\rho\\alpha\\beta = O.'t 

10 addSymmetry (' W_\\mu\\nu\\alpha\\beta', 1, 0, 2, 3) 

11 def w = ' W''\\mu\\nu_\\alpha'^\\beta = W^\\mu\\nu_\\alpha'^\\bat . t 

13 def N = 'N'^WrhoWalphaWbeta = O'.t 

15 def M = 'M_\\alpha'^\\beta = M_\\alpha^\\beta.'t 

17 def F = 'F_\\mu\\nu\\alpha\\beta = F_\\mu\\nu\\alpha\\betat 

19 def div = oneloopdiv4 (Kmv , K, S, W, N, M, F) 

21 def counterterms = Elim±nateFromSymmetr±es>> div . counterterms 

23 counterterms = ' M'* \ \mu_\ \mu = M'.t >> counterterms 

24 counterterms = ' W_\\mu\\nu''\\alpha_\\alpha = W_\\mu\\nu.'t >> 

counterterms 

25 counterterms = ' W^\\alpha_\\alpha = W. t >> counterterms 

27 println counterterms 

[> counterterms = -M+ ( 2 / 3 ) *F_{ \nu\beta } { \epsilon }_{ \rho_{ 5 } } *F^ { \ 
nu\beta\rho_{ 5} }_{\epsilon} - (32/135) *R_{ \mu\nu } *R'^ { \mu\nu 
}+ (44/135) *R**2+ (1/9) *R*W- (1/9) *R_{ \mu\nu} { \mu\nu } + ( 1 /2 4 ) 
*W'^ { \ epsilon\delta\alpha}_{ \ rho_{ 5 } } *W_{ \epsilon\delta}''{\ 
rho_{5} }_{\alpha}+ (1/48) *W_{ \delta } " { \delta\alpha }_{ \rho_ 
{5} }*W_{\beta}"{\beta\rho_{5} }_{\alpha} 



3 
4 
5 
6 



+ 'd_\\alpha'^\\beta!rl/3*(' + 
' g'^\\mu\\nu«rg''\\gainina\\delta ' + 
'+ g''\\mu\\gamma»:g''\\nu\\delta ' + 
'+ g'^WmuWdeltatg'^WnuWgamma) ') .t 
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Multiplying this result by l/167r((i — 4) and integrating over the space-time volume gives: 



~ " 167r(d - 4) y ^ V 3 ' ' 



135 135 9 9 ^ 



where the notation M^'^ = M, Wf^^^a = and W"a = W used. As in the case of the second 
order operator, it is necessary to explicitly specify, that tensors S and N from (6) are equal to 
zero. 



5 Basic internal architecture 

5.1 Mappings of indices 

Perhaps the most significant difference between tensor- and symbol-oriented computer algebra 
systems lies in the comparison of mathematical expressions. In the symbol-oriented CASs the 
result of atomic comparison problem^ is just a logical true or false, while in the case of tensor- 
oriented CAS it transforms into a complex pattern matching problem, which produces a complex 
object as a result. The comparison problem of tensorial expressions will be revealed in this 
section through several examples. 

The most common question, which can be asked about two expressions, is whether they 
define the same tensor (to within a free indices relabelling). This question arises in such frequent 
routines like substitutions and reduction of similar terms. Consider the following expressions: 

FabG'' Fi.G'^^ 

These two expressions have the same tensorial structure. If rename a to i and c to j in the left 
expression one will get exactly the same tensor as defined by the right expression. So, the result 
of such comparison is not just true or false, but a mapping of free indices of one expression onto 
free indices of another expression. 

One of the major quirks of the problem lies in the fact that free and dummy indices hold 
completely different places. Mappings of free indices are global for expression, while dummy 
indices have their scopes. If some free index is present in several places, then its mapping will 
be the same everywhere. Consider the following two sums: 

FabG\. + MadN'', FigC^, + M.giV?, 

Here a and c should be renamed into i and j respectively in both summands in order to transform 
l.h.s. into the r.h.s. If such renaming does not exist, then two expressions define different tensors: 

FabG'c + MadN\. -h FigG'^j + MjgN", 

On the other hand dummy indices have more sophisticated mapping rules, since their explicit 
names are not important (consider indices b and d in the l.h.s and q in the r.h.s.). So, dummy 
indices have their scopes, which introduces additional complexity into the problem. 

The formal definition of mapping of indices can be formulated as follows: the mapping of 
indices of tensor A to tensor B is a reflection between free indices of tensor A and indices of 

^determination of whether two expressions are equal, i.e. operation that is the main building block of such 
complex routines as pattern matching 
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tensor B (need not be free), such that if we rename indices of A according to it, then A becomes 
exactly the same tensor as B (to within dummy indices relabehing). 

In fact, mapping of indices is the most frequently executed "atomic" operation in any real 
calculation. The implementation of mapping procedure and its underlying abstraction is a key 
achievement of Redberry. 



5.1.1 Multiple mappings and symmetries of tensors 

In general, several mappings of indices can exist for a pair of tensors. Consider the following 
primitive example: 

Zij + Zjj > Zab + -^fea 

As one can see, there are two possible mappings: 

It is clear that mapping of tensor onto itself gives permutational symmetries of its indices. So in 
case of the above primitive example, one can find that: 

gives 

where Mi corresponds to identity symmetry, while M.2 means that if T^j = Zij + Zji, then 
T — T 



5.1.2 Sign of mapping 

Mappings have a sign property: 

Zij — Zji — y Zab — Zba 

gives 

fz— )'ol , ( i ^ b 

Mi = + < . , > and 7W2 = - i . 

Negative sign of mapping A^2 means that in order to obtain the r.h.s., one needs to apply 
mapping to the l.h.s. and negate it. This is an important addition to the definition given in the 
beginning of this section. 

Consider a more complex example. Suppose that tensor Rab is antisymmetric, then: 

gives 

Ml = + { b-y j ) and M2 





Sign property of mappings and processing both symmetries and anti-symmetries in a common 
way makes these entities in Redberry fully consistent with each other. 
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5.1.3 Mappings and substitutions 



Usually mappings are used inside a complex routines. In most cases mapping rules are calcu- 
lated to be applied to some tensors. Common pattern of mappings usage can be illustrated by 
their application in substitutions. Suppose one need to apply the following substitution rule: 

to the expression: 

In order to do this, one need to perform the following steps: 

51. Find the l.h.s. of the substitution in the expression: Fp^G^'^ 

52. Perform the comparison, i.e. build mapping of indices: 



a ^ p 



S3. Check for dummy indices clash and resolve it. Here tensor R^ap^i. conflicts with FprG^'^, 
so we can resolve problem by adding additional mapping rule: 




54. Apply mapping to the r.h.s. of the substitution: 

55. Replace: 

FprG^ y R praG^ 

5.1.4 Examples and performance 

Comparison of mathematical expressions is the most frequent atomic operation in the program 
and its performance directly affects execution time of any calculation. Because of this, its 
implementation is one of the main targets for optimization. 

The first remarkable thing concerning performance, is that in some cases (e.g. when tensors 
have a lot of symmetries etc.) the number of all possible mappings can be huge (about A^!, 
where is a number of free indices). However, in most cases there is no need to compute 
all possible mappings. Thus, Redberry uses iterator pattern to enumerate possible mappings, 
and the calculation of each subsequent mapping occurs only on corresponding step of iteration. 
Consider example: 

1 setAntiSymmetric (' R_ah' ) 

2 def from = 'R_{ab}*A_c + R_{bc}*A_a' .t, 

3 to = 'R_{ij}*A_k + R_{ jk}*A_i' .t 

4 def mappings = from % to 

5 mappings . each { println it } 

[> +{_a -> _i, _b -> _j, _c -> _k} 
[> -{_a -> _k, _b -> _j, _c -> _i} 
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Here, mappings object is an iterator over possible mappings from tensor from to tensor to , 
which calculates mappings only on demand. It has the main method take () (++ in Groovy 
syntax) that calculates and returns the next mapping or null if no more mappings exists. 
Consider the following examples: 



1 def from = 'A_m^n'.t, to = 'A^a_b'.t 

2 (from % to) . each { 

3 println it 

4 assert to == (it >> from) } 
> +{_m -> -^a, '"n -> _b} 

5 addSymmetryC R_ahcd' , 0, 2, 1, 3) 

6 from = ' R_a''bi_l ' . t 

7 to = 'R_aj'^b'' j'.t 

8 (from % to) .each { 

9 println it 

10 assert to == (it >> from) } 

[> +{_a -> _a, -> '~b, ''i -> _j, _1 -> ^j} 
[> +{_a -> _a, ^b -> _j, ''i -> ^b, _1 -> '"j} 

11 setAntiSymmetric (' A_mn' ) 

12 setAntiSymmetric {' F_itinah' ) 

13 from = ' (A_m^n - A_m^p«rA_p^n) *F_nk^i_ j + A_mn»:A^n_ j*A''i_k' . t 

14 to = '-(A_d''a + A_p^a(tA_d'^p)*F'^d_kq''i - A^a_bkA'^b_q*A''i_k' . t 

15 (from % to) .each { 

16 println it 

17 assert to == (it >> from) } 

[> +{"1 -> _k, _j -> _q, _k -> ^i, _m -> '"a} 
[> -{''i -> '^i, _j -> _q, _k -> _k, _m -> '~a} 



As it seen, it is possible to map indices with different states, i.e. upper indices onto lower and 
vice- versa (for metric types only). In order to apply mapping rules to tensor one can use >> 
operator (as in lines 4, 10, 18). This will automatically perform raising or lowering of indices if 
it is meant by mapping and resolve dummy clashes. 

5.2 Graph theory and tensors 

As was mentioned in the previous section, finding of possible mappings between indices of 
tensors is a computationally hard problem. One of the main ideas, used in the mappings module 
in Redberry, is that contractions of indices define a mathematical graph. 

For simplicity, let's assume anywhere where it is not specified, that all tensors in this section 
are completely symmetric. Here is a formal description of representation of a product as a 
mathematical graph: 

• Each multiplier matches a vertex of the graph which label is equal to tensor name (see 
below) 

• Each pair of dummy indices matches an edge of the graph 

^It implements an OutputPortUninterruptible pattern from rcdberry-pipe library for concurrent 
algorithms. 



37 



• Each free index matches a special vertex, labeled ' — ', and connected to the corresponding 
multiplier. So, free indices are in some sense contracted with abstract multipliers. 

Here, term "tensor name" is used to distinguish the mathematical nature of different tensors. 
In case of simple tensors two tensors are considered to have different names if and only if they 
have different mathematical nature. For example, all of the following tensors 

A A ^ A 

are considered to have equal names in Redberry, but, for example, names of tensors Bp^y or 
A^^^ differs from those listed. Internally, Redberry represents simple tensor names as integer 
numbers and assigns unique integers for simple tensors with different mathematical nature. If 
some multiplier is a complex tensor (e.g. sum), then Redberry uses its hash code as a label for 
graph vertex. 

In order to illustrate proposed rules, let's consider the following tensor: 

This tensor has two free indices a and (3, and since all of the multipliers are completely symmetric, 
it is easy to show that this tensor is also symmetric: 

Aa'^'^ Afj^fj^Cpy = Ap^'^ Apa^C py. (13) 

Fig. (la) and (lb) illustrate graphs that represent the l.h.s. and r.h.s. of the above equality. Since 
these tensors are equal, they have the same graph representation, which is shown in Fig. (Ic). 




(c) Universal graph representation 



Figure 1: Graph representation of tensors Aa^'^Ap^^Cpy and Aj^^'^ Apa^Cpy. Figures (a) and (b) 
show correspondence between graph edges and indices of tensors. Figure (c) shows the common 
graph representation. 

It is clear that symmetry (13) leads to a nontrivial automorphism of the corresponding graph 
which is shown in Fig. 2. On the other hand, nontrivial graph automorphisms, which permute 
free indices (' — ' vertices), are equivalent to tensor symmetries. 
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Figure 2: Automorphism of the tensor graph A^^^ A^p^Cpy, which proofs its symmetry proper- 
ties. 

The above example gives a good idea of the tight connection between the tensor-oriented 
computer algebra and computational graph theory. Many problems specific to tensor-oriented 
CAS can be re-formulated or reduced to the known problems in the graph theory and vice-versa. 
It is now clear that mappings problem is closely related with graph isomorphism problem. 

In order to take into account particular permutational symmetries of simple tensors, Redberry 
performs some king of an additional edge labeling. So, in general case of arbitrary tensor, we 
have disconnected, undirected graph with labeled vertices and edges. Graph problems, such as 
effective finding of isomorphisms or automorphisms of graphs are one of the central problems in 
computer science, and actively studied for decades (see [14]). Many algorithms were proposed to 
solve these problems, and while many of them seem to perform well on random graphs, a major 
drawback of these algorithms is their exponential time performance in worst cases. Moreover, it 
is still an open question: whether graph isomorphism (or automorphism) is V or AA'P-complete. 

Additionally to the graph isomorphism problem, more complex problems arise when building 
mappings between tensors, which maps free indices onto dummy indices. Consider the following 
example. Suppose we need to build a mapping from tensor F to tensor T, where 



Tensor F have two free indices ^ and A, while tensor T is a scalar. One can see, that mapping 
should simply contract free indices of like 



Nevertheless, this simple answer involves building a special kind of mapping between the cor- 
responding graphs, which are shown in Fig. 3. We see, that this mapping is not an exact 
isomorphism, but a surjection. It is clear, that calculation of such mappings implies finding 
of isomorphic structures in corresponding tensors, and can be effectively reduced to the graph 
isomorphism problem. 

Another computationally hard problem arises when applying substitution rules with a product 
in the l.h.s. Suppose one have a substitution rule A ^ where j4 is a product of tensors, and 
one needs to apply it to tensor C . Then one needs to find all subproducts of C, such that there 
is mapping from it to A. This implies, that one needs to find all subgraphs of graph C which 
are isomorphic to graph A. This is a well-known AAP-complete subgraph isomorphism problem 
[15], and algorithms known to date, in general, have an exponential or even factorial execution 
time in the worst case. 



from tensor F: A^o^^App" As^f^ A^'^^ A'^'i'^ A^^^^A^'"'' A''. 

i 

to tensor T: A.^.^App't' A^^^ A^'^^A'^^'' A^^^A^^"" AP 



At 
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Figure 3: A complex example of mapping construction. Here the left graph represents tensor 
"from" with free indices ^ and A. The right graph represents tensor "to", which is scalar. The red 
dashed lines shows the mapping between isomorphic structures, while black dashed lines shows 
the mapping from "from" free indices onto target (dummy) indices in "to". 

Redberry has internal implementation of graph and subgraph isomorphism problems. The 
implementation is based on the topological hashing and hash refinement mixed with depth-first 
search. It actively utilizes additional information about symmetries of tensors, and the worst 
case occurs only when there are a lot of similar symmetric multipliers in the target products. 
While, the algorithms developed in Redberry have a pretty good performance in the most cases 
arises in physical problems, they still have bad performance in some cases, so improvement of 
these algorithms is one of the main goals of current development of Redberry. 

5.3 Expression-tree traversal and modification 

As mentioned in Sec. 2.2 each tensor in Redberry is a container of its child tensors, so any 
complex expression becomes a hierarchical tree of such containers. Iteration over direct child 
elements of a tensor described in Sec. 2.2. Besides, there are a special tools for iteration and 
modification over a whole tree. 

There is a core class ^ in Redberry that performs traversal over any given expression. It 
basically generates a sequence of traversal events like: entering or leaving of subexpressions. 
The main feature of this class is its ability to in-place modify a tree while traversing. Redberry 
provides simple facade classes for tree traversal. There are two basic ways to perform a depth-first 
search on expression tree: 



1 def t = 'a + Sin[x + y] '.t 

2 t . parentAf terChild { a -> print a.toStringO + ', ' } 

[> a, y, X, y+x, Sin[y+x], a+Sin[y+x], 

^ TreeTraverselterator in package cc . redberry . core . tensor . iterator 
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3 t . parentBef oreChild { a -> print a.toStringO + ', ' } 
[> a+Sin[y+x], a, Sin[y+x], y+x, y, x. 



This example illustrates steps of parent-after-child and parent-before-child iteration modes. It 
is also possible to specify a particular guide in order to restrict iteration only to a parts of 
expression^ : 



4 t = 'a*b + c + Sin[x + y]'.t 

5 def guide = { 

6 tensor, parent, positionlnParent -> 

7 if (tensor. class == Sin) ShowButNot Enter 

8 else if (parent == 'a*b'.t) DontShow 

9 else Enter 

10 } as TraverseGuide 

12 t . parentAf terChild (guide) { a -> print a.toString() + ', ' } 

[> c, b*a, Sin[x+y], c+b*a+Sin [x+y ] , 

13 println( ) ; 

14 t . parentBef oreChild (guide) { a -> print a.toStringO + ', ' } 



[> c+b*a+Sin [x+y ] , c, b*a. Sin [x+y]. 



The guide takes three arguments: current tensor, its parent tensor and integer position in parent 
tensor and returns TraversePermission . In such a way guide tells the iterator to step over 
current tensor ( DontShow ), or show current tensor, but not to go deeper into its children 
( ShowButNotEnter ), or to continue normal iteration ( Enter ). The guide shown above guides 
the iteration as follows: if current tensor is of type Sin , then it will be shown during iteration, 
but the iterator will not go inside it; if parent of the current tensor is equal to a*b , then such 
tensor will not be shown at all; otherwise iterator will continue normal iteration. 

Modification of expression tree can be performed in the same way. As an example, lets 
consider a naive implementation of a substitution transformation: 



1 def subs = { 

2 t, expr -> 

3 t . transf ormParentAfterChild { 

4 def c = ++(expr[0] % it) 

5 if (c != null) c >> expr[l] 

6 else it 

7 } 

8 } 

10 def t = ' z_m*Cos [x_m*y''m - x_n*(z''n + fn) ] + t_m'.t 

11 println subs(t, ' z_a + t_a = y_a'.t) 



[> y_in 



^one should import 


Sin 


from 


cc . redberr .core. tensor. functions 


7 


TraverseGuide 


from 


CO . redberr . core . 


iterator 


and also statically import 


TraversePermission 


from the last pack- 



age. 
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The key aspect of parent-after-child modification, is that if child node is changed, then its parent 
node will be reduced to standard form (see Sec. 2.4) before it will be shown^. 

Another important thing concerning expression tree modification, is that any conflicting 
dummy indices will be relabeled automatically: 



12 t = 't_a*(x + t_a*t''a*x) ' .t 

13 println subs(t, 'x = F_a''a'.t) 

> t_{a}* (F_{b}"{b} + t_{b}*t"{b}*F_{c}"{c}) 



6 Conclusions 

In this paper we presented Redberry — an open source computer algebra system designed 
to manipulate with symbolic tensorial expressions. Redberry is a computer algebra system, 
which considers both tensors and indexless expressions in a common way. It provides basic 
tensor-specific routines such as tensor symmetries, multiple index types, dummy indices handling, 
M^H^X-style I/O, mappings of tensor indices and a set of tensor-specific transformations etc. 
Rich functionality of Redberry was demonstrated on a complex physical problems: Feynman 
graphs and one-loop counterterms calculations. Compared to many other tensor-oriented CASs, 
Redberry provides a simple and convenient facade to modern high-level programming language 
(Groovy), which makes it possible to use all features of general-purpose programming language 
combined with domain-specific features required by computer algebra. 

At the moment, Redberry does not implement many features, which should be implemented 
in a general-purpose tensor-oriented CAS. However, the architecture of Redberry is built to be 
open for further extensive development. This is a first public release of Redberry CAS, and the 
development is ongoing. We currently working under the following issues: pattern matching, 
derivatives, expressions standard form improvement, performance of symmetries, mappings, and 
factorization. 

Redberry code is well covered by over 800 unit tests. However, there are several known 
bugs (that not affects correctness of calculations), which are listed in Redberry issue tracker. 
Redberry performance can be demonstrated by the following numbers^: the whole calculation 
from Sec. 4.1 (Compton scattering in QED) takes about 0.5 seconds, calculation of counterterms 
of non- minimal squared vector field takes less then 3 minutes, non-minimal vector field — about 
22 seconds, minimal fourth order operator — about 2 seconds. The calculation of trace of 8 
gamma matrices takes less then 20 ms, while, for example, FeynCalc makes it in 300 ms (on 
the same hardware). The calculation of trace of 12 gamma matrices takes less then 3 seconds in 
Redberry. 

Redberry is licensed under GNU GPLv3. Additional documentation, examples and installa- 
tion instructions are available at http://redberry.cc. The source code and issue tracker can be 
found at http://bitbucket.org/redberry. We are welcome for contributions and any suggestions. 

^ The first modification will be performed on the sum z^n + t '"n , which will be replaced with y'^n . After 
this modification, the whole branch of the expression tree, which contains this term, will be permanently reduced to 
the standard form. So, when iteration will reach the argument of cosine, which will be x_m* y^m - x_n*y'^n , 
it will be reduced to zero. On the next step the cosine receives zero as an argument and reduces to 1. On the last 
step of the iteration, it finally becomes possible to replace obtained z_m + t_m with y_m . 

^on Intel Core 15 2.27 GHz 
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A Inputting matrix expressions 

Matrices are tensors with a subset of indices marked as matrix indices. These indices have 
nonmetric types, which means, that raising and lowering are forbidden. The rules of matrix 
multiplication can be formulated as follows. Suppose we have two tensors 



I and 5^"i-"fc"r r 



where Ui and /j denotes matrix upper and matrix lower indices respectively, and /i and denotes 
a whole subset of other indices. Then the product Tf^^ = A^By has the following form: 

= Am"'-"""c,...c„,5/^"""'""'+i-"'=",,..4^ if ni < K, 
y «i...«„„ ^ 4 «i...«n„ r if ni>ku, 

f^'^ Iku + l-^nih-hi Cl...Cfc„tfc„ + l...(„; U h...lki ' 

These definitions becomes especially clear in a particular cases. Consider, for example, three 
matrices: matrix with one upper matrix index (usually called vector), matrix 7^*^- with one 
upper and one lower index, and matrix Vi with one lower matrix index (usually called covector). 
Then using the above rules, one can obtain the well-known expressions: 

T^u = 'l^llu T^u'j = Ifi'du^'j (matrix) 

= IfiV T^u' = Ifi'cV'' (vector) 

= ^^7m = Vcl^^i (covector) 

T^u = V j^-fu V T^u = Vi-f^'c-fu^jV^ (scalar) 

Redberry has a tool to specify what tensors should be considered as matrices. It automatically 
inserts the additional matrix indices according to the formulated rules. Consider the examples: 

1 defineMat rices 'G_a', 'G', Matrixl. matrix, 

2 'v', Matrixl. vector, 'cv', Matrixl. covector 

4 println 'G_a*G_b'.t 

[> G_a"a'_b' *G_b'^b'_c' 

5 println 'G_a*v'.t 

\> G_a"a'_b' ★v'^b' 

6 println 'cv*G_a'.t 

[> cv_a' *G_a'"a'_b' 

7 println ' G*G_a = G*v*cv*G_a + f_a'.t 

> G'^a' c'*G a'^c' b'= G'^a' c' *v^c' *cv d' *G a'^d' b' +f a*d''a' b' 
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8 println ' cv*G_a*G_b*v + g_ab'.t 

[> v_a' *G_a''a'_b' *v^b' + g_ab 

9 println 'Tr[G_a*G_b + n_b*G_a] + n_a*n_b'.t 

t> G_a'"a'_b' *G_b'"b'_a' + n_b*G_a'"a' _a' + n_a*n_b 
10 println 'Tr[G_a*G_b + n_b*G_a + n_a*n_b]'.t 

[> G_a^a'_b' *G_b''b'_a' + n_b*G_a^a' _a' + n_a*n_b*d^a' _a' 



In the first two lines we specify tensors which will be considered as matrices in further calcu- 
lations. Method def ineMatrices takes a set of string representations of matrices, where the 
corresponding matrix signature (type and number of upper/lower matrix indices) is specified 
after each set. By default, Redberry provides four predefined nonmetric types: Matrixl (Latin 



lower case letters with strokes). 


Matrix2 


(Latin upper case letters with strokes). 


Matrix3 


(Greek lower case letters with strokes), Matrix4 
number of upper and lower indices of matrix type 


(Greek upper case letters with strokes). The 
is specified using the corresponding property: 


vector 


means one upper index. 


covector 


means one lower index. 


matrix 


means one upper 



and one lower index. In the general case of p upper and q lower indices, one can do 



11 def ineMatrices 'M_\\mu', Matrixl. tensor (2, 3) 
13 println ' G*M_\\alpha' . t 

[> G"{a' }_{f' }*M_{\alpha}"{f'b' }_{c'd'e' } 



As we see from the lines 7 and 10, Redberry inserts similar free matrix indices for the l.h.s. 
and r.h.s. of the expression. Moreover, if, for example, some sum contains both matrix (with 
non empty free matrix indices) and non-matrix expressions (without matrix indices), then the 
later will be multiplied by the identity matrices, i.e. Kronecker deltas (compare to line 8). 

Redberry provides a special syntax for traces of matrices shown in the lines 9 and 10. By 
default the trace is taken with respect to all matrix indices, but if some expression is a matrix 
with respect to several indices types, then it is possible to specify particular types of indices for 
which the trace should be taken: 



1 def ineMatrices 'A, B', Matrixl. matrix, Matrix2. matrix 

3 println ' Tr [A*B] ' .t 

[> A-^ia' }_{b' }"{A' }_{B' }*B'^{b' }_{a' }"{B' }_{A' } 

4 println 'Tr[A*B, Matrixl] '.t 

[> A"{a' }_{b' }"{A' }_{C' }*B"{b' }_{a' }"{C' }_{B' } 



B Redberry artifacts structure 

Redberry project consists of three main artifacts, which are published in Maven Central 
Repository: 
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• cc . redberry. core — the core framework written in pure Java, which contains the 
implementation of basic Redberry routines and general-purpose transformations. 

• cc . redberry . physics — applications of Redberry in real physical problems. It contains 
tools for Feynman graphs and one-loop counterterms calculation. 

• cc . redberry. groovy — a Groovy facade for Redberry. It contains the classes and cate- 
gories^ defining DSL features and syntax notations for common routines from redebrry-core 
and redberry-physics . 

The current version of all Redberry artifacts is 1.1. All official releases of Redberry artifacts 
including binary jar files, sources, documentation and tests, can be found in Maven repository. 
Redberry sources and version control are hosted at bitbucket.org/redberry. Redberry is free and 
open source project and anyone can contribute to Redberry with new features using Bitbucket 
forks. 



^see Groovy documentation 
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